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Abstract. The collective ground-state correlations stemming from low-lying quadrupole excitations are 
computed microscopically. To that end, the self-consistent mean-field model is employed on the basis 
of the Skyrme-Hartre-Fock (SHF) functional augmented by BCS pairing. The microscopic-macroscopic 
mapping is achieved by quadrupole-constrained mean-field calculations which are processed further in 
the generator-coordinate method (GCM) at the level of the Gaussian overlap approximation (GOA). We 
study the correlation effects on energy, charge radii, and surface thickness for a great variety of semi-magic 
nuclei. A key issue is to work out the influence of variations of the SHF functional. We find that collective 
ground-state correlations (GSC) are robust under change of nuclear bulk properties (e.g., effective mass, 
symmetry energy) or of spin-orbit coupling. Some dependence on the pairing strength is observed. This, 
however, does not change the general conclusion that collective GSC obey a general pattern and that their 
magnitudes are rather independent of the actual SHF parameters. 
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1 Introduction 

The nucleus is a highly correlated many-body system and 
thus the topic of correlations has always accompanied nu- 
clear physics. Short-range correlations dominate due to 
' the huge short-range repulsion in the nuclcon-nucleon in- 
teraction. Their treatment requires very involved approa- 
ches as, e.g., Brueckner-Hartree-Fock, variational Jastrow, 
hyper-netted chain, or no-core shell model calculations, 
see e.g. [TJEIEEEI ■ Large-scale applications employ sim- 
pler approaches, among which the most microscopic ones 
are self-consistent mean-field models based on effective 
energy-density functionals or effective forces, respectively. 
The functional is motivated by a formal derivation from an 
underlying many-body approach, see e.g. [6], but adjusted 
phenomenologically because of the as yet insufficient in- 
put from ab-initio models; for a recent review see [7]. A 
widely used functional is, e.g., the Skyrme-Hartree-Fock 
(SHF) approach. However obtained, the nuclear energy- 
density functional is supposed to contain all crucial cor- 
relations effectively. The reasoning is quite similar to that 
of density-functional theory in electronic systems where 
theorems and clean theoretical developments have come 
much farther, for an overview see [8]. The perfect func- 
tional thus allows describing a system by a mere mean- 
field calculation and yet obtaining the correct energy and 
density. That ideal functional is known to be highly non- 
analytic to account for the sudden changes from one par- 
ticle number to another [8]. Actual functionals, however, 



are rather smooth functions of the involved densities, par- 
tially owing to approximations in their derivation and par- 
tially for practical reasons of manageability. This, in turn, 
means that the functional can incorporate only the smooth 
trends of the correlations and will probably fail to ac- 
count for quickly changing contributions. Yet such effects 
are known to exist in the nuclear landscape. The nuclear 
shell structure leads to dramatic changes in the nuclear 
shape when adding or removing nucleons, known as the 
nuclear Jahn- Teller effect (9lll0j. Pairing moderates these 
strong fluctuations [9|lllj , but sizeable fluctuations remain 
as can be seen from the changes in the low-energy spec- 
tra along isotopic or isotonic chains [12]. The lowest ex- 
citation in even-even nuclei is in most cases a quadrupole 
state assigned as 2^. This 2^ state is a collective state. 
It is associated with substantial rccoupling of the simple 
mean-field excitations which, in turn, produces significant 
correlation effects. These collective correlations are not 
smooth with changing particle number and thus are not 
included in the energy-density functional. They have to 
be added explicitly. That has been discussed for decades, 
but mostly based on semi-empirical estimates, see e.g. [Ml 
HP!] . Thorough calculations based on self-consistent mean- 
field models have also been tried for a long time; see e.g. 
[T5 l fT6 l fr7 l H8j. Extensive surveys running over all nuclear 
landscapes based on SHF have been published recently 
19 , 20, 21j and there exists also a systematic study using 
the Gogny force [22] showing very similar trends. It is the 
aim of this paper to give a survey of collective correla- 
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tions within the SHF mean field, exploring in detail the 
dependence on the SHF functional. The aim is to prepare 
ground for a revision of the phenomenological adjustment 
of SHF by finding out the least correlated nuclei. 

The low-lying quadrupole states can be formulated in 
terms of the Bohr-Hamiltonian, which establishes collec- 
tive dynamics in the five quadrupole degrees of freedom 
[2"5] . The parameters of the collective Hamiltonian are 
usually adjusted phenomenologically, see e.g. the appli- 
cations in [24| and more recently in the interacting boson 
model |25j . The direct connection between a microscopic 
description and a collective picture can be established by 
a collective path, which consists of a continuous series of 
mean-field states with prescribed deformations, produced 
by a quadrupole constraint. The mapping of the path into 
collective dynamics is established by virtue of the gen- 
erator coordinate method (GCM). [26ll27j . The complex 
integral equations of the GCM can be simplified using the 
Gaussian-Overlap- Approximation (GOA) which, further- 
more, allows to establish contact between the microscopic 
foundation and the collective Bohr-Hamiltonian [281IT7]. 
In practice, there are, on the one hand, fully fledged GCM 
calculations which skip the collective Hamiltonian as in- 
termediate level and compute low-energy spectra directly 
from the coherent superposition of the collective path; 
these sophisticated calculations imply exact projection for 
the conserved quantities, like particle number, angular 
momentum, and center of mass; from the many published 
results we mention here [291I30] as two recent examples. 
On the other hand, there are the techniques which use the 
Bohr-Hamiltonian as an intermediate stage, for an early 
example see |15pi6| and for more recent achievements [3 1 1 
152] . We are going here to use the GCM-GOA. The path is 
generated from axially deformed SHF-BCS states. The full 
five-dimensional topology is regenerated by interpolation 
into the triaxial plane in an extension of the method pre- 
sented in [33j . The interpolation scheme promises reliable 
results for vibrators whose collective deformation rarely 
reaches large values. It is thus well suited for the present 
survey dealing with chains of semi-magic nuclei. 



2 Formal framework 

2.1 The Skyrme mean-field model 

The starting point for the self-consistent microscopic de- 
scription is the SHF energy functional E = E(p,T,J,j,a), 
which is expressed in terms of a few local densities and 
currents obtained as sums over single-particle wave func- 
tions: density p, kinetic density t, spin-orbit density J , 
current j, and spin density er, each occurring twice, once 
for protons and once for neutrons. It is augmented by a 
pairing functional deduced from a zero-range two-body 
force (DI=delta interaction) or combined with a density 
dependence (DDDI=density-dependcnt delta interaction) 
|34j . They are described with the pairing functional 
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where DDDI-pairing corresponds to po = 0.159 fm~ 3 , while 
Dl-pairing is recovered for pa — > oo. The pair density 
is x( r ) = J2 a u a v a\fa{ r )\ 2 ■ The real part thereof con- 
fines the contributions to time-even parts. This coincides 
with the standard pairing functional for stationary states 
(where x is purely real) and it implies that time-odd pair- 
ing effects in dynamical response are ignored. Pairing is 
treated within the BCS approximation and augmented by 
the Lipkin-Nogami (LN) correction which serves to force a 
minimal amount of pairing everywhere and so prevents the 
sudden changes associated with the pairing phase transi- 
tion. That smoothing is compulsory for the further pro- 
cessing of the mean-field states in the description of col- 
lective dynamics. The BCS approximation is applicable 
for well bound nuclei where the Fermi energies are safely 
below the continuum threshold. This holds for the nu- 
clei considered here along the whole collective deformation 
path. For details about SHF, BCS and LN, see the review 
article [7J. 

The mean-field equations are derived variationally from 
the given energy functional. They determine the ground 
state or a locally stable isomer. Collective motion goes 
through a succession of deformed mean fields, the collec- 
tive path. The low-lying quadrupole mode is particularly 
soft and the collective path is related to quadrupole de- 
formed shapes. We generate the path by imprinting a ded- 
icated deformation using quadrupolc-constrained mean- 
field equations, often called constrained Hartree-Fock (CHF), 



^MF = huF — — XQ20 
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where /imf is the SHF-BCS mean field Hamiltonian (de- 
pending on the local densities and thus on the state |^ Q20 )), 
£f is the Fermi energy, and A the Lagrange parameter for 
the quadrupole constraint. Q20 is the quadrupole operator 
and a 2m its dimensionlcss expectation value, i.e. 
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with A the total particle number and r the r.m.s. radius. 
The index m can run over —2, — 1, 0, 1 and 2. For a while, 
we will consider only m = which corresponds to axially 
symmetric deformations. The states are labeled with the 
dimensionlcss quadrupole moment (|2d[) . The quadrupole 
operator is modified by a Wood-Saxon-like damping func- 
tion / cu t with an extension of three times the nuclear ra- 
dius and width of 1 fm in order to suppress the unbound 
regions from the asymptotics oc —x 2 ,—y 2 , or — z 2 |35j . 
The equations are solved with an extra iterative loop to 
maintain a wanted value of a 2 o |36j . This is done for a 
dense set of deformations a 2 o- At the end we have the 
collective path {|^a 2 o)} as a series of mean- field states 
along which the collective motion can evolve. Its further 
processing will be discussed in section 12.21 Before carry- 
ing on, we ought to mention that the optimal scheme for 
generating a collective path is adiabatic time-dependent 
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Hartrcc-Fock (ATDHF) which generates the constraints 
in a self-consistent manner without imposing a precon- 
ceived idea of the collective deformation, see e.g. ]28[[37] . 
Experience shows that the quadrupolc constraint is a good 
approximation in cases where the low lying 2^ state is a 
truly collective mode. That applies practically to all nu- 
clei, except for the doubly-magic ones. Thus our results 
near doubly-magic nuclei have to be taken with care. 

The SHF functional sets only a framework within which 
there is a manifold of different parameterizations. Most 
available parameter sets describe ground states properties 
equally well but differ in other observables like, e.g., exci- 
tation spectra or nuclear matter properties [7]. One thus 
should consider several sufficiently different parameteriza- 
tions to distinguish particular features of a given parame- 
terization from general features of SHF. We will consider 
the following standard Skyrmc parameterizations: SkM* 
as a widely used traditional standard [38] , Sly6 as a recent 
fit which includes information on isotopic trends and neu- 
tron matter [39], and SkI3 as a recent fit which maps the 
relativistic iso- vector structure of the spin-orbit force [40] . 
The set contains distinct effective masses where SkM* has 
m* jm = 0.8 while SLy6 and SkI3 have significantly lower 
values (0.69, or 0.6 respectively). SkI3 differs in that it 
has basically no proton-neutron coupling in the spin-orbit 
mean field. Thus all three forces differ somewhat in the 
actual shell structures they produce. Besides the effective 
mass, the bulk parameters (equilibrium energy and den- 
sity, incompressibility, symmetry energy) are comparable. 
In addition, we employ a set with a systematic variation of 
effective mass, symmetry energy and spin-orbit force |41j 
to explore the sensitivity to these aspects of the force. 



2.2 The collective Hamiltonian 

Large-amplitude collective motion proceeds along the col- 
lective path, a series of mean-field states which differ very 
little in energy. Its description requires a coherent super- 
position of the states along the collective path. That is 
done in the framework of the GCM. The GOA allows 
to map the microscopic picture into a collective Hamil- 
tonian for the five quadrupole degrees of freedom a = 
(«2-2, Q!2-i! Q!20! Q!2i, 022) [28]. The final result is summa- 
rized briefly. 

The collective Schrodinger equation to be solved reads 
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with the collective Hamiltonian 



H c ° n -- i {V M ,{B^(a),V,}} + y(a) . (4) 

It is obtained from the general rules of a collective map 
which apply analogously to the collective map iV co11 of the 
particle- number operator [28] . The details are given here 
for the Hamiltonian. The potential in H co11 is given by the 
mean-field energy from which the spurious contributions 



from rotational, vibrational, and center-of-mass zero-point 
energies are subtracted, i.e. 
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where the are the generators of a change in the quadru- 
pole momentum. The inverse collective mass tensor is eval- 
uated as in terms of the constraint mean-field Hamiltonian 
Hmf and the 
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where Q M are the generators of a dynamical boost. They 
are determined sclf-consistcntly as dynamical linear re- 
sponse to the collective momentum V v , i.e. 
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where Hyyp is the mean-field Hamiltonian belonging to 
state {'Pa) and h lcS p,Q is the response Hamiltonian for 
time-odd perturbations about |# a ) [42]. The mass thus 
obtained is called self-consistent or ATDHF mass. It is 
to be noted that most straightforward GCM calculations 
omit the dynamical linear response and deal with an ap- 
proximate mass which can be deduced purely from the 
static configurations. That amounts to replacing the dy- 
namical boost operator in formula ([6a| simply by the 
redundant complement to [28], i.e. 



(7) 



Accounting for the dynamical response explores the dy- 
namic path for slow collective motion. Within the assump- 
tions of the GOA framework it thus corresponds to a 
dynamic generator-coordinate method (DGCM). We will 
compare later on these two approaches to the collective 
mass. The collective path |<P a ) is computed for axially 
symmetric shapes, i.e. along a^o- This determines the raw 
collective potential -Eshf, the masses B 00 and the widths 
Aoo along the a^o- By a slight rotation of the nucleus, 
generated by the angular momentum operator J x or J yi 
the moments of inertia are evaluated providing additional 
information about rotational properties of the nucleus. In- 
cluding this information allows to evaluate the and 
Ai,-! components of the GOA tensors. The data from pro- 
late (ck20 > 0) and oblate {a.20 < 0) configurations arc in- 
terpolated into the triaxial plane with methods similar to 
those employed in [33] but now accounting for the topol- 
ogy of the five dimensional quadrupolc configuration space 
(see appendix [XJ. The requirement of rotational symme- 
try for the collective Hamiltonian H coU and of smoothness 
in the triaxial direction allows to determine all needed in- 
gredients of the collective Hamiltonian in the full triaxial 
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plane, for the regime of small deformations. It is to be 
noted that the assumption of small deformations limits 
the present treatment to nearly spherical nuclei. These 
are the ones which wc will study in the following. 

Coming back to the collective Schrodingcr equation 
([3]), we sec that it contains a constraint on particle num- 
ber. Each mean-field state |^a 20 ) is tuned to the same 
average particle number (N) = N by tuning the Fermi 
energy in the mean-field Hamiltonian (|2bj) . But the coher- 
ent superposition may lead to a small drift of the average 
particle number. The term oc <5ep in ((3]) serves for a final 
particle- number correction (PNC). 



3 Results and discussion 

3.1 Ground-state correlations of semi-magic nuclei 

FigureCD shows the results of the SHF plus DGCM-GOA 
calculations for the correlation effects on the four bulk 
observables energy E, charge r.m.s. radius r, diffraction 
radius R, and surface thickness a for all semi-magic nuclei 
considered in this survey. The upper panel complements 
that by the energies of the low lying 2 + compared with 
the experimental results (taken from [13]). That figure is 
generic in the sense that the trends seen here are typical 
for a great variety of Skyrme forces. It summarizes the 
correlation effects from low lying 2 + states. 

The uppermost panels show the excitation ener- 
gies. The heavier systems distinguish nicely the doubly 
magic nuclei by large E(2f ). The experimental energies 
drop suddenly when going away from the doubly magic 
stage while the theoretical results show a somewhat softer 
transition. That qualitative mismatch is cause by two ap- 
proximations which are not optimal for doubly-magic nu- 
clei: the LN recipe for stabilized pairing and a simple 
quadrupole constraint rather than full ATDHF. It is to 
be noted that the same effect was found in the study of 
[2Tj which differs in several details of the treatment but 
does also employ a collective path from CHF. There are 
slight quantitative discrepancies in the mid-shell regions. 
These change with the Skyrme force and can be related 
to the level density at the Fermi surface. Such variations 
will be discussed later. 

The correlation effects (middle and lower panel) show 
the expected trends, being minimal at doubly magic points 
and growing large at mid-shell. The question of interest is 
how large they actually become. The horizontal lines indi- 
cate typical values of acceptable uncertainties in a well fit- 
ted Skyrme force, I MeV for the energy and about 0.02 fm 
for charge radii or surface thickness. The lighter systems 
generally have larger correlation effects, often beyond the 
desired bounds. Some Z=20 and N=20 systems even have 
slightly positive correlation energies, with 40 Ca being the 
worst case. That unphysical effect is due to the GOA. It 
requires smooth collective paths [44]. But for these small 
semi-magic systems, the abrupt pairing phase transition 
lies well inside the range of the shape fluctuations of the 
ground state and it is not sufficiently well smoothed by 
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Fig. 2. Collectivity of the low lying 2 + states for Sn isotopes 
(left) and N = 82 isotones (right) computed with SkI3 and 
DI pairing. Upper panels: comparison of E(2f) excitation en- 
ergies, from pure mean field (two quasiparticle energies), full 
collective calculations, and experiment. Middle panels: B(E2)] 
values from collective calculations and experiment. Lower pan- 
els: correlation energies. 

the LN scheme. An E COIT +0.5 MeV is anyway within 
the expected precision of the method. The 40 Ca is more 
dramatic because the low 2 + state is not really collective. 
Besides these critical light nuclei, correlation effects are 
negative for energies and positive for r.m.s. radii, as it 
should be. The parameters derived from the charge form 
factor, diffraction radius and surface thickness can go both 
ways and do so. The results carry two surprises: first, the 
form parameters (middle panels) generally show smaller 
correlation effects relative to the goal than the energies, 
and second, the correlation energies arc large within iso- 
topic chains while remaining moderately small in isotonic 
chains. Both features call for a revision of fitting strate- 
gics. A thorough discussion of reliable data sets for the 
adjustment of Skyrme parameters will follow in a separate 
publication. Here we continue with discussing variations 
of Skyrme forces, pairing, and collective approximations 
to corroborate the above result. To that end, we confine 
the presentation to the isotopic Z=50 chain and the iso- 
tonic N=82 chain. That contains all necessary trends while 
rendering the figures more transparent. 

3.2 Low lying 2+-excitations 

Figure [2] shows more details about the 2+ state and the 
correlations energies. The upper panel contains the E(2f) 
excitation energies and adds the two-quasiparticle (2QP) 
energies E^qp for comparison. The difference between the 
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chains of semi-magic isotopes chains of semi-magic isotones 

Z=20 Z=28 Z=50 Z=82 N=20 N=28 N=50 N=82 N=126 




2028 28 4050 66 82 110 126 142020 28 40 5050 58 66 82 90 
Neutrons N Protons Z 

Fig. 1. Correlation effects on semi- magic isotopic (left) and isotonic (right) chains computed with Skl3 and DI pairing. Upper 
panels: excitation energies, theoretical and experimental as indicated. Middle panels: correlation shifts of charge r.m.s. radii 
r, diffraction radii R, and surface thicknesses a. Lower panels: correlation energies. 



Eiqp and the fully dressed E(2f) characterizes the strength 
of the residual interaction which is clearly related to the 
correlation (lowest panels). One also sees the systematic 
effect that isotones (upper right) have less collective down- 
shift than isotopes (upper left), and thus smaller correla- 
tion energy (lowest panels). 

The middle panel complements the excitation energies 
by the corresponding B(E2) values. These agree in the av- 
erage fairly well with the experimental data (taken from 
g3llii71iHllini[^[ni[511S31). However, there are substan- 
tial deviations in details: a dramatic drop at N=66 and 
68 and a trend to overestimation towards shell closure 
N=82. The breakdown of the B(E2) in the mid-shell tin 
isotopes is related to the kink in the E(2f ) which, in turn, 
can be traced back to a kink in the two-quasiparticle en- 
ergy (upper panel). The reason is the interference with a 
shape transition which is probably unrealistic, an artefact 
of that particular force (we see it also for SLy6, but not for 
SkM* and the other forces discussed here). The differences 
around the doubly-magic tin isotope result from using the 
CHF path as an approximation to the optimized ATDHF 
path. By construction CHF overestimates the collectivity 
of the ground state and thus also the transition probabil- 
ity to the 2 + -state. The mismatch is similarly seen in the 
too large excitation energies. The same artefact appears 



also in studies based on different nucleon-nucleon interac- 
tions using a similar collective model based on a CHF path 
[2"T] while a recent QRPA study performs much better in 
that respect [M] . As QRPA can be considered as the small 
amplitude limit of ATDHF, that result indicates that full 
ATDHF will cure the mismatch for those isotopes. 

The lowest panels show the cm. correlation energies 
and pairing energies in addition to the collective corre- 
lation energy. The cm. energy has a very smooth trend 
because it is related to excitation across shells. The col- 
lective correlation energies, on the other hand, are related 
to intra-shcll excitations and thus changing from isotope 
to isotope. The trend is very similar to the trend of the 
pairing energy which depends similarly on the available 
intra-shell phase space. 

Figurc[3]dcmonstrates the effect of the collective ground 
state correlations (GSC) on the trend of the deviation 
from experimental binding energies. The pure mean-field 
calculations (open circles) display an unresolved trend for 
the isotopic chain (left) while having already the correct 
trend for the isotones (right). The GSC now perfectly com- 
pensate that due to strong correction for the mid-shell iso- 
topes (see figure^]). The practically flat trend of the GSC 
for the isotones leaves their already agreeable pattern in- 
tact. 
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Fig. 3. Difference between theoretically computed and ex- 
perimental binding energies in Z=50 isotopes (left) and N=82 
isotones (right) computed with SLy6 and DI pairing. Com- 
pared are pure spherical mean-field calculations (open circles) 
with correlated ground state energies (filled circles). 
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Fig. 4. Effect of the particle-number correction on excitation 
energies (top) and on correlation energies (bottom) for Sn iso- 
topes (left) and N = 82 isotones (right) computed with Skl3 
and Dl pairing. 



3.3 Comparing approximations 

Figure [4] shows the effect of particle number correction 
(PNC) on excitation and correlation energies. The case 
without PNC means to ignore the term oc Sep in The 
effect on E(2^) is small, but the PNC are slightly improv- 
ing the steepness of the transition from doubly-magic nu- 
clei to semi-magic ones. The effect on correlation energies 
is still small but noteworthy at a quantitative level. It is 
particularly interesting that it always reduces E COTI while 
the general trends remain unchanged. The PNC reacts 
particularly sensitive at N=66 and 68, the region where 
shape fluctuations lead to the much reduced B(E2) val- 
ues. 

Figure [5] deals with the treatment of the collective 
mass, comparing results of the (standard) GCM which 
employs only stationary deformation paths yielding the 
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Fig. 5. Effect of the dynamical response ("dynamic GCM") on 
excitation energies (top) and on correlation energies (bottom) 
for Sn isotopes (left) and N = 82 isotones (right) computed 
with SkI3 and DI pairing. 



redundant mass with the dynamic GCM which ex- 
plores the dynamical response of the nucleus and so em- 
ploys the self-consistent cranking mass (|6bj) . The dynamic 
GCM lowers the E(2+) and improves, in particular, the 
steepness for the step away from doubly-magic nuclei. The 
B(E2) values react particularly sensitive at N=66 and 68, 
precisely the region which is very sensitive for some forces. 

For the correlation energies, the trends remain robust 
qualitatively, but there is a sizeable reduction by about 
20-50%. Both trends can be explained by the feature that 
the dynamical response is a first step towards the fully 
self-consistent ATDHF path (that is why self-consistent 
cranking is often called ATDHF cranking). The variation 
in the larger space reduces the excitation energies and, at 
the same time, takes a different cut through the collective 
landscape which eventually reduces correlations [55] . This 
shows the importance of proper dynamical response and 
indicates that full ATDHF will improve the performance 
at and around doubly magic nuclei. 



3.4 Variation of the Skyrme- and pairing force 

Figure [6] compares results for three widely used Skyrme 
parameterizations. The global trends of the results are 
the same: smaller E COTT for isotonic chains, overall size of 
Eqott, growth of E corT towards mid shell, and too smooth 
decrease away from doubly-magic nuclei. There are differ- 
ences in detail. The dip in the E(2^) and B(E2) at 116 Sn 
does not appear for SkM* , and this force also has the cor- 
relation maximum at a different position in the isotopic 
chain. 
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Fig. 6. Correlation energies and 2 + -energies calculated with 
SkI3, SLy6 and SkM* and Dl-pairing. 



It is not possible to judge which of the features of the 
three forces is most responsible for the slight differences. In 
order to disentangle the driving agents, we now are consid- 
ering a series of Skyrme parameterizations produced un- 
der the same conditions while varying systematically one 
selected feature. It originated from a systematic analysis 
of giant resonances and their dependence on the Skyrme 
parameters [H] . Figure [7] shows a variation of the effec- 
tive mass m*/m. This feature has an impact on the shell 
structure, mainly the level density. The E(2~£) are almost 
inert. The B(E2) values are also generally inert, but show 
some sensitivity around the notoriously critical N=66. The 
same holds for -E C orr where some effects can be seen in re- 
gions around N=60 and N=82 for Sn isotopes and here the 
correlations increase with decreasing m* /m as one would 
expect. The other regions where E co „ is less sensitive are 
dominated by pairing as we will see later. The variation 
of the symmetry energy shown in figure [8] is very robust 
and shows even smaller effects. 

We also studied other variations of bulk parameters 
(sum rule enhancement k, incompressibility, density de- 
pendence of a sym ). Changes in those cases were even smaller 

Figure [9] shows the effect of using different forms for 
the 1-s coupling [41]. Recall that there is some variety in 
handling the 1-s forces for SHF functional. The standard 
form is derived from a zero-range two-body 1-s interaction 
[56] which leads to a fixed mix of isoscalar and isovec- 
tor 1-s coupling, the isovector being half as large as the 
isoscalar one. That is indicated by "stdls" in the figure. 
The classical limit of the relativistic mean- field model sug- 
gests that the isovector spin-orbit should be nearly zero 
|40j . a variant which is indicated by "rmfls" in the plot. 
Moreover, there is a spin-orbit term oc J 2 emerging from 
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» a m*/m = 1,oq 

50 55 60 65 70 75 80 85 50 55 60 65 70 
neutron number N proton number Z 

Fig. 7. Influence of the effective nucleon mass m*/m on ex- 
citation energies (top) and on correlation energies (bottom) 
for Sn isotopes (left) and N = 82 isotones (right). Compared 
are results for three Skyrme forces with systematically varied 
m* /m while keeping all other features the same. 



Z=50 isotopes N=82 isotones 




50 55 60 65 70 75 80 85 50 55 60 65 70 
neutron number N proton number Z 



Fig. 8. Influence of the symmetry energy a sym on excitation 
energies (top) and on correlation energies (bottom) for Sn iso- 
topes (left) and N — 82 isotones (right). Compared are results 
for four Skyrme forces with systematically varied a sym while 
keeping all other features the same. 
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Fig. 9. Influence of the form of the l*s functional on excita- 
tion energies (top) and on correlation energies (bottom) for Sn 
isotopes (left) and N = 82 isotones (right). 




50 55 60 65 70 75 80 
neutron number N 



85 50 55 60 65 70 
proton number Z 



Fig. 10. Influence of the pairing strength and pairing func- 
tional on excitation energies (top) and on correlation energies 
(bottom) for Sn isotopes (left) and N = 82 isotones (right). 
Compared are results for DI pairing with systematically varied 
strength and DDDI pairing, all for SkI3. 



the kinetic zero-range interaction. It is called tensor 1-s 
term. Many parameterizations omit that term. Inclusion 
is indicated in the figure by "tens" and omission by "no 
tens" . a sym The impact of the different 1-s models is well 
visible: Inclusion of tensor 1-s yields somewhat larger cor- 
relation energies and lower E(2 + ) at places, particularly 
in combination with the standard 1-s model. The B(E2), 
on the other hand, are robust towards shell closures and 
show some sensitivity mid shell. It is obvious that the 
1-s force has a large effect on shell structure. We observe 
that changing the 1-s model can even change the ordering 
of levels near the Fermi energy. This, in turn, modifies 
quasi-particle energies and thus E{2 + ) as well as E co „, 



Figure [TOl demonstrates the effect of pairing. It shows 
a variation of the strengths V$ q of the pairing functional 
([U for the case of Dl-pairing relative to the original DI 
and DDDI pairing functional. The effects are considerable 
for all three observables shown. Again, the E(2f) are very 
sensitive due to a close breakdown of pairing for the lower 
strength. It is noteworthy that DDDI even produces differ- 
ent trends than DI for E(2f). The correlation energies are 
more robust, yet showing here the largest variations of all 
cases studied here. The results for the E(2^) in compari- 
son to experiment show that the reduced pairing strength 
is clearly ruled out. What remains proves again that the 
trends of correlation energies and their approximate mag- 
nitudes are a very robust quantity, a generic feature of all 
reasonable SHF functionals. 
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Fig. 11. Comparison of different approximations to the 
quadrupole correlated state. The DGCM-GOA and GCM- 
GOA calculations made use of the topological interpolation 
to recover the full quadrupole configuration space (5D) during 
the solution of the GCM-GOA Hamiltonian. The GCM cal- 
culations are performed directly by coherent superposition of 
the microscopic states along the static axial path (3D) avoid- 
ing the intermediate step of the GOA approximation and the 
collective Schrodinger equation. 



3.5 Comparison with a large-scale GCM analysis 

Figure [IT] tries a comparison with the GCM calculations 
of [I9jl20j. It has to be taken with care because there are 
several differences at once: The method of [T51I2T)] handles 
pairing in a slightly different phase space (cutting symmet- 
rically above and below the Fermi energy), it deals with 
full particle number projection (where we restore particle 
number in the average), it uses full angular momentum 
projection (while we do that in GOA), it uses only the 
stationary path, it performs the collective superposition 
only within the space of axially symmetric shapes (where 
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Fig. 12. Correlation diagram between 15 corr and the pairing 
energy iJ pa i r of the static mean-field ground-state for isotonic 
chains (upper panel) and isotopic chains (lower panel) . The cor- 
relation and pairing energies were calculated using the Skyrme 
parameterizations SkI3, SLy6 and SkM* and Dl-pairing. The 
average linear trend E COII = 0.25-E pa i r is indicated as the solid 
line through the scatter while typical deviations of ±lMeV are 
indicated by the dashed lines. 



we interpolate to fully triaxial collective dynamics), and 
it circumvents the collective Schrodinger equation by de- 
termining the superposition weights directly through the 
Hill- Wheeler equation with the given overlap kernels. Fi- 
nally, we deal with two slightly different forces, SLy4 and 
SLy6. But these forces belong to the same family and be- 
have very similar such that this should not spoil the com- 
parison. In spite of all differences, the results in figure [TT1 
arc quite instructive. First of all, there is an overall agree- 
ment in trends and magnitude. Looking more closely, one 
sees even a close agreement between the full, but static, 
GCM treatment of [l9j[20] and our approach at the same 
level, namely the one using the mere GCM mass (denoted 
"GCM" in the plot). The main difference seems to come 
from the self-consistent mass in DGCM. That statement, 
however, has to be taken with care because the multitude 
of differences in treatment leaves many options open. It 
would require a very detailed study to disentangle the var- 
ious ingredients. This is beyond the scope of the present 
paper. 



3.6 A ro 



ugh estimate of the average correlation energy 

The question arises whether one can establish a relation 
between the correlation energy and other quantities which 
are simpler to compute. Pairing properties are the most 
promising candidates for that purpose because pairing 
is switched on and off with similar trends as collective 
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Fig. 13. The reduced correlation E rc d = E COII — 0.25-E pa i r 
drawn versus the proton number Z for the semi-magic isotone 
chains (upper panel) or neutron number N for the semi-magic 
isotones (lower panel). As in figure [12] the data were obtained 
with the three Skyrme parameterizations SkI3, SLy6 and SkM* 
and DI pairing. 



correlations. We have figured out that the closest con- 
nection exists with the pairing energy ([I]) of the mean- 
field ground state. Figure [12] collects results from a va- 
riety of forces in a diagram £? C orr versus -Epair- The up- 
per panel shows isotopic chains (i.e. neutron pairing) and 
the lower isotonic ones (i.e. proton pairing ) as indicated. 
The nuclei with large pairing and correlation energies in- 
dicate a clear linear trend which we have determined as 
-Ecorr = 0.25 Ep^ by least squares fitting. Small nuclei 
(open symbols) show a larger scatter, but have smaller cor- 
relation energies throughout. The Ni isotopes with Z=28 
and the isotones with N=28 are particularly off the line. 
The fit line is accompanied by two parallel lines indicat- 
ing the band ±1 MeV. It is gratifying to see that most 
nuclei with N, Z > 50 stay fairly well within that band. 
That nourishes the hope that one may develop a simple 
formula for previewing correlation effects. It ought to be 
mentioned, however, that the trend is not so strictly linear 
for DDDI pairing. This, as well as the strong deviations 
for the N=126 isotones here, have yet to be understood. 

Taking the simple trend from figure 1121 we define a 
reduced correlation energy as i^od = E corT — 0.25 -Epair. 
Figurc[T3]shows the trends of -E rc( j over all isotopic and iso- 
tonic chains. That shows even better the nice performance 
for the heavier nuclei and some yet unresolved trends. The 
latter call yet for refinement of the estimate and, at the 
same time, may serve to give a clue for the improvement. 



4 Conclusions and outlook 

Based on nuclear mean-field calculations using the Skyrme 
energy-density functional, we have studied ground state 
correlations as produced by the soft modes associated with 
low-lying quadrupole states. The collective dynamics was 
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described by a coherent superposition of deformed mean- 
field states which were produced by quadrupole constrained 
Skyrme-Hartree-Fock calculations. The information on col- 
lective transport (masses, inertia) was evaluated by a self- 
consistent calculation of the dynamical response to chang- 
ing deformation. The thus correlated state was handled 
with the generator-coordinate method (GCM) in the Gaus- 
sian overlap approximation (GOA) yielding effectively a 
collective Hamiltonian in quadrupole coordinates and sim- 
ilar mapping of other obscrvablcs. The study concentrated 
on semi-magic nuclei with more or less shape fluctua- 
tions about the spherical shape. The microscopic mean- 
field calculations employed axially symmetric deforma- 
tions while the collective dynamics was considered with 
all five quadrupole degrees of freedom by properly inter- 
polating the collective properties (potential, masses and 
quantum corrections) into the triaxial plane. The proce- 
dure was employed for a systematic exploration of ground 
state correlations on nuclear bulk properties (energy, radii, 
surface thickness) scanning its dependence on the underly- 
ing energy functional and on detailed corrections as, e.g., 
restoration of particle number or dynamical response. We 
have also investigated the energies and B(E2) values of 
the low-lying 2 + excitation. 

On the side of the model development, we find that 
the proper treatment of dynamical response has moder- 
ate effect on E(2^) excitation energies (about 10-20%) 
but a huge one on correlations (factor two on E COIY ). The 
restoration of the correct particle number makes less ef- 
fects, a few percent on E(2^) and about 10% on E COII . 
We also find that the low lying quadrupole state is a 
truly large-amplitude mode. Even in doubly magic 208 Pb, 
a small amplitude limit turns out to be unreliable (15% 
effect on E(2f)). Some small nuclei cause problems in 
that they have a positive correlation energy. The GOA 
becomes insufficient there due to a small number of nu- 
clcons involved in the motion (lack of collectivity) and 
due to a strong impact of the pairing phase transition. 
For system sizes, the trends of E(2^) do not reproduce 
the experimentally found steep decrease when going away 
from doubly-magic nuclei. This defect was also found in 
earlier calculations using GCM. That is a problem of the 
simple quadrupole constraint used for generating the col- 
lective path. More elaborate schemes as, e.g., adiabatic 
time-dependent Hartrcc-Fock are required to describe cor- 
rectly the transition from collectivity in mid-shell nuclei 
to predominantly two-quasiparticle excitations in doubly- 
magic systems. 

The large scale variation of energy functionals has shown 
that the correlation effects are very robust. Their trends 
and magnitudes come out quite similarly for all Skyrme 
forces in this survey. Even the more sensitive 2 + excita- 
tion shows very similar patterns throughout. Differences 
are seen, of course, in detail. Variation of bulk proper- 
ties as symmetry energy or effective mass has little effect. 
Somewhat more sensitivity is seen for variations of the 
spin-orbit model, in particular tensor spin-orbit can pro- 
duce up to 20% on E COII and 40% on E(2f). The strongest 
effect comes from variation of the pairing strength. Par- 



ticularly the E(2^) are very sensitive which is no surprise 
because of known strong interplay between shape transi- 
tions and pairing. However, the basic result of the robust- 
ness of -E corr remains valid even here. We have also found 
a close relation between E COII and the pairing energy of 
the mean-field ground state which show in the average a 
linear dependence on each other. 

The results provide clear indications for the next steps. 
On the formal side, one needs to go for the self-consistent 
constraints according to ATDHF replacing the simple- 
minded quadrupole constraint and one should search for a 
better smoothing of the pairing phase transition, if possi- 
ble. For further applications, one can take advantage of the 
robust information on correlation effects and use these as 
guide for searching least correlated ground state observ- 
ables as input for better adjustments of Skyrme energy 
functionals. It is also worth while to improve the esti- 
mate of correlations energies in terms of simple to compute 
quantities. Work in all these directions is in progress. 
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A Interpolating the Collective Hamiltonian 

A.l Intrinsic system 

The Cartesian representation of the quadrupole coordi- 
nates (a M ) is advantageous for formulating the GCM-GOA 
Hamiltonian For the solution of the collective Schro- 
dinger equation a formulation in terms of the intrinsic 
coordinates (/3,7,#) is more convenient [23 1 I5T 1 I24] : 



i9{cos(7)D^(*)+sin(7)- 



>(#) 



V2 J • 

(8) 

The shape of the nucleus is parameterized in one of the 
principal axes systems of the nucleus by the total deforma- 
tion f3 and triaxiality 7 while the orientation of this prin- 
cipal axis system is given by the three Euler angles 1?. The 
representation of the collective Hamiltonian in these coor- 
dinates is advantageous due to the decoupling of rotation 
and vibration modes. In detail the collective Hamiltonian 
reads 

6= V t B(/3, 7 )V + y(/3, 7 ) , (9a) 
with collective potential V, mass tensor 
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and the gradient 

iL' x y iX^ 

- I Op, ^a 7 , 2 /3sin( 7 -^) ' 2/3sin( 7 -^) ' 2/3sin( 7 ) 

(9c) 

It is important to note that the components of the col- 
lective mass tensor as well as the potential depend on 
the deformation coordinates 7) but not on the ori- 
entation of the intrinsic system ( r d). The actual orien- 
tation of the system is incorporated by the angular mo- 
mentum operators of the intrinsic system L' k = L' k {'&) 
(k G (x, y, z) = (1,2,3)). The choice of the intrinsic sys- 
tem is not unique, which results in additional symmetry 
conditions: 



(ft 7) 



3) 



k 2n 



(10a) 
(10b) 



with X e {V, Bpp, B 77 , S/3 7 } which have to be fulfilled 
necessarily by the collective parameter functions [57] . 



obtained from the axial calculations: 
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where Jf( p ) = X(a 20 = /?) and = A(a 20 = -/?). By 
inserting the expansion functions i?i and Vi into (| 1 2 j) the 
collective mass tensor and potential arc found. 



A. 2 Expansion of Collective Potential and Masses 



A. 4 Zero-Point Energies 



Following similar lines to the rotationally invariant expan- 
sion of the generalized Bohr-Hamiltonian in [57] a sym- 
metry conforming, polynomial expansion of the collective 
potential and mass tensor can be derived by appropriate 
coupling of the spherical tensors of coordinates a and con- 
jugate momenta 7r. 

To account for the fluctuations of the collective poten- 
tial with deformation j3 while assuming a smooth inter- 
polation in the triaxial degree of freedom, the polynomial 
expansion of the collective potential V is approximated 
up to lowest order in 7 while maintaining symmetry re- 
quirements and the most general /3-dependency similar to 
[33]: 

708,7) = Vi{fi) +/3 3 cos(37)Vi(/?) . (11) 

By means of similar arguments the 7-smooth parameteri- 
zation of the collective mass tensor is given in ([T2]) . To end 
up with a polynomial expansion the functions Vi (/3) , Bi (/3) 
have to be analytic and even functions. 



A. 3 Interpolation Formulas 



The zero-point energies (ZPE) are composed from the in- 
terpolated potential V{(3, 7), the mass tensor B(/3, 7) and 
the GOA width tensor A(/3, 7). Thereby the width tensor 
is interpolated in the same manner as the collective mass 
tensor. In the-five dimensional case the zero-point energy 
correction also includes terms from spurious ZPE in the 
additional degrees of freedom, namely the 7- vibration and 
^-rotation mode as well as contributions from the (3-j- 
coupling. The contributions from the kinetic ZPE arc con- 
structed straight forward from the parameter functions by 



77 



kin 



ZPE 1 \ jdiiv 



(14) 



^ppBpp + \ 11 B 11 + 2\p~ ( B{3~ ( +y^Afci3fc 
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A compact notation of the potential ZPE requires to de- 
fine the curvature tensor of the collective potential = 
^VjuVj/V with intrinsic components 



The parameterization of the collective potential ([TT]) and 
mass ([12]) are adjusted to the values obtained by CHF and 
cranking along the axially symmetric cuts (a = a 2 o, a ±2 = 
a ±1 = 0). Besides the Born-Oppenheimer surface V(ct2o) 
the three dimensional axial GCM-GOA Hamiltonian is de- 
fined from the (3- vibrational mass -8/3/3(0120) corresponding 
to motion along the axial symmetric path and the degener- 
ate rotational mass B x / y (a2o) related to rotations around 
the non-symmetry axes. Generalizing the collective Hamil- 
tonian to the five dimensional configuration space yields 
the simple interpolation rules between the expansion co- 
efficients in (fTS"]) and the collective parameter functions 



Cp 1 

c k 



fa ~ ^ 



V(j3,i) 



co t( 7 -fcf) 

— p — *> 



^,7) ■ (15) 



In terms of the components of the curvature tensor C the 
potential ZPE is evaluated in analogy to the kinetic ZPE 



12 



P. Kliipfel et al: Systematics of collective correlation energies from self-consistent mean-field calculations 



%(^,7) = Bi{p) +/3cos(3 7 )B 2 (/3) +P 2 B 3 {f3) + /3 3 cos(3 7 )B 4 (/?) , 

BpyiP.l) = -PM*Y)Ba(P) " ^-sin(3 7 )B4(/3) , 

B 77 (/3, 7 ) = Bi(/3) -/3cos(3 7 )B 2 (/3) + ^B 3 (/9) + y cos(3 7 )B 4 (/3) , 

B,(/3,7) = -/3cos( 7 )B 2 (/3) - ^ (7 - 10cos(2 7 )) S 3 (/3) - £ (7cos(3 7 ) - 10cos( 7 )) B 4 (f3) . (12) 



and reads 
1 



E 



ZPE 



pot — 2 ^ ^ 



(16) 



k=x,y,z 



The topological interpolation scheme (fTTl [T2")) enforces the 
asymptotic behavior for the ZPEs when approaching the 
spherical limit: As expected from the harmonic approxi- 
mation the contributions from the two vibrational modes 
and the three rotational modes degenerate yielding five 
equivalent contributions for (3 = 0. This feature of our 
model enables an improved description of ground state 
properties of magic and especially doubly-magic nuclei. 
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